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Abstract 

The thermodynamic, kinetic and magnetic properties of the hydrogen monomer on doped graphene layers 
were studied by ab initio simulations. Electron doping was found to heighten the diffusion potential barrier, 
while hole doping lowers it. However, both kinds of dopings heighten the desorption potential barrier. 
The underlying mechanism was revealed by investigating the effect of doping on the bond strength of 
graphene and on the electron transfer and the coulomb interaction between the hydrogen monomer and 
graphene. The kinetic properties of H and D monomers on doped graphene layers during both the annealing 
process (annealing time to =300 s) and the constant-rate heating process (heating rate or =1.0 K/s) were 
simulated. Both electron and hole dopings were found to generally increase the desorption temperatures 
of hydrogen monomers. Electron doping was found to prevent the diffusion of hydrogen monomers, while 
the hole doping enhances their diffusion. Macroscopic diffusion of hydrogen monomers on graphene can 
be achieved when the doping-hole density reaches 5.0 x 10 13 cm . The magnetic moment and exchange 
splitting were found to be reduced by both electron and hole dopings, which was explained by a simple 
exchange model. The study in this report can further enhance the understanding of the interaction between 
hydrogen and graphene and is expected to be helpful in the design of hydrogenated-graphene-based devices. 

PACS numbers: 68.65.Pq, 67.63. -r, 68.43.Bc 
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I. INTRODUCTION 



In recent years, there have been many investigations into hydrogenated graphite surfaces 
and graphene due to their importance in astronomic exploration, 1 ^ nuclear industry, 4 - 5 graphite 
(graphene)-based hydrogen storage 6 and graphene-based electronic devices 

Study of the thermodynamic and kinetic properties of hydrogen adatoms on graphene can help 
to understand the interaction mechanism between the hydrogen adatom and the graphite surface 
(graphene layer). The usability of the newly proposed hydrogenated-graphene-based devices relies 
heavily on these properties. Many experiments have been carried out on the thermodynamic and 
kinetic properties of hydrogen adatoms on graphite surfaces , 2 ' 3 ' 19 ' 20 substrate-supported graphene 
layers (monolayer and multilayer) 8 ' 9 ' 11 ' 13 ' 14 and even the free standing graphene layer. 9 21 The hy- 
drogenation of graphene starts with the adsorption of a hydrogen monomer, which breaks an aro- 
matic n bond in graphene and makes the adsorption of other hydrogen atoms very easy. Thus, the 
thermodynamic and kinetic properties of the hydrogen monomer are essential to the hydrogenation 
process, which are having been investigated theoretically in recent years.^2 2 .- 8 The ab initio simu- 
lation from our previous paper 28 has given some predictions of the kinetic properties of the hydro- 
gen monomer on neutral graphene, which closely reproduced some experimental observations^. 
Furthermore, the hydrogen monomer has been used by some theorists to engineer the electronic 
structures, magnetic properties and transport properties of some graphene-based devices.- 3 ^ 1 The 
thermodynamic and kinetic properties of the hydrogen monomer are critical to the stability of 
these devices. 

The charge doping of graphene is a common phenomenon in intercalated graphite 
compounds, 32 34 metallic surface-supported graphene layers^ 5 - - — and gated graphene layers.^ 3 ^ 
Charge (electron or hole) doping was expected to influence the thermodynamic and kinetic prop- 
erties of hydrogen adatoms (including monomers) on graphene. In addition, the magnetism in 
graphene-based materials has recently drawn tremendous interest from the scientific community 
due to their lower density compared to transition metals, compatibility with biological systems, 
plasticity and so on, 40 which introduces another set of applications for charge-doping effects on the 
hydrogenated graphene. Overall, charge doping provides a possible approach to modifying var- 
ious properties of hydrogenated-graphene-based materials and devices to fulfill the requirements 
for various applications. 

In this report, the thermodynamic and kinetic properties of the hydrogen monomer on doped 
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graphene layers were simulated using a composite method consisting of density functional the- 
ory (DFT), 41 density functional perturbation theory (DFPT) 42 and harmonic transition state theory 
(hTST). 43 " 45 The mechanism of the charge-doping effect on the thermodynamic and kinetic prop- 
erties of the hydrogen monomer was revealed by investigating the electronic structures and phonon 
spectra. The kinetic properties of H and D monomers on doped graphene layers during both the 
annealing process (to =300 s) and the constant-rate heating (1.0 K/s) were simulated. The effect of 
charge doping on the magnetic properties derived by DFT calculations have been well explained 
by a simple exchange model. 

H. METHODOLOGY 

The adsorption energy of the hydrogen monomer E ads is defined as the energy difference be- 
tween the totally desorbed state and the adsorbed state 

E ads = E H + E GL - E GL+H (1) 

where E H , E GL and E GL+H are the total energies of an isolated hydrogen atom, an isolated graphene 
layer and a graphene layer with a hydrogen monomer adsorbed on it, respectively. Electron doping 
will make the anti-bonding n* orbital in graphene occupied by electrons, while hole doping will 
reduce the electronic occupation of the bonding n orbital. Thus, both the bonding n and the anti- 
bonding n* orbitals, which are together described by a total n bond, are considered here. The 
chemisorption of a hydrogen monomer will break an aromatic II bond. The breaking energy of a 
IT bond (En) is related to the adsorption energy of the hydrogen monomer in the chemisorption 
state, according to Ferro's analysis. 48 Therefore, a simple model is used to reexpress E ads as 

E ads (cr) = E C - H ((T) - E n (cr) (2) 

where <x is the doping-charge density and E c _ H (cr) is the formation energy of a C-H bond. If 
E C - H is assumed to be unchanged under charge doping, the contribution of E n (cr) to E ads (cr) is 
explicitly addressed in this way. The adsorption energy based on this assumption is expressed as 

E* ads (o-)=E C - H (0)-E n (o-) (3) 

Comparing E ads in Equ. © with E*. in Equ. © gives the relative contributions of E c _ H (cr) and 
E n (cr) to E ads (cr). E u can be calculated from the adsorption energies of the hydrogen monomer 
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and para-dimer 48 



Eu(cr) = E p 2 a {cr) - 2E ads (cr) (4) 

where E pa J s a is the adsorption energy of a hydrogen para-dimer on graphene. 

The over-barrier jump frequency v between two local minimum states (initial and final states or 
reactant and product states) is calculated from the quantum-mechanically modified hTST , 28 ' 43 ' 46 ' 47 
in the Arrhenius form of 

F 

where v* m is the exponential prefactor and E ac is the activation energy. The activation energy E ac 
is defined as the vibrational zero-point energy corrected potential barrier, which is expressed as 

3JV-1 -, 3N 



Eac = Ay p + i£^f-i£^ 

i=\ i=l 

= AV P - AF vib (0) (6) 



where AV P is the potential barrier in the reaction path, which can be obtained from DFT calcula- 
tions; cjl and a)? are the vibrational frequencies of the ith mode in the initial and saddle-point (SP) 
states, respectively, which can be obtained from DFPT calculations; AF vih (0) is the vibrational 
zero-point energy correction. The total vibrational degrees of freedom are 3N. An imaginary vi- 
brational mode along the migration coordinate in the SP state is excluded from calculation. Thus, 
there are 3N-1 vibrational modes considered for the SP state. The quantum-mechanically modified 
prefactor is expressed as 28 



3N ~ l luo s 



k B T ! 



I! exp(T^)n r (o;f) 



=i 



Vqm ~ M,- , i\ 

i=i 

where h T ((x>i) is the bosonic phonon occupation number of the ith vibrational mode. 

The first-order rate equation for the desorption of the hydrogen monomer is defined %g&2M± 

dn(t) 

—j— = -v des (T)n(t) (8) 
at 

where Vd es is the desorption jump frequency of the hydrogen monomer; n(t) is the residual number 
of hydrogen monomers on a graphene layer at time t\ t = represents the starting time of the 
kinetic movement. 



4 



In the annealing process (T = T ), the variation of the residual number with respect to the 
annealing time interval (t ) is expressed as 

n(t) = n(0)exp[-v des (T )t ] (9) 

where n(0) is the monomer number at the starting time (t = 0). Conceptually, the relative residual 
monomer number (n(t)/n(0)) is also the desorption probability of a single monomer. In this case, 
the diffusional property can be described by the mean square displacement (< |r(f) - r(0)| 2 >) of a 
monomer parallel to the graphene layer. From Fick's second law, we have 

<\r(t)-r(0)\ 2 >=2dD ad (T )t (10) 

where d is the dimensionality of the diffusion of a hydrogen monomer on graphene, taken to be 2 
here, and D ad is the temperature-dependent diffusion coefficient of the adatom, described by^ 

D ad = ^-Ya 2 (11) 
2d 

where Y is the total jump frequency of the monomer, and a is the jump length. For the diffusion 
of a hydrogen monomer on graphene, Y is taken as 3v di ff (y d iff is the diffusion jump frequency) 
because there are three equivalent paths for the diffusion of a monomer on graphene, and a is taken 
as the optimized C-C bond length of 1 .426 A. The diffusion radius of the monomer is defined as 
the square root of the mean square displacement 



r dif = yjldDJMt (12) 

which directly determines the diffusional mobility of the hydrogen monomer on graphene. 

In the constant-rate heating process (T = at), the variation of the residual number with respect 
to time is 



n(t) =n(0)exp[- f v des (T)dt] (13) 
Jo 

where the heating rate a is always taken as 1.0 K/s in experiments . 2i 19 ' 20 In this case, the diffusion 
radius is expressed as 



r dif = \\2d D ad (T)dt (14) 



In this study, a hydrogen monomer on a 5 x 5 periodic supercell of graphene (50 C atoms), with 
a 10 A vacuum along the direction perpendicular to the surface, was taken as an isolated monomer. 



The DFT and DFPT calculations were carried out using the Quantum Espresso code package. 49 
The ultrasoft 50 spin-polarized PBE 51 pseudopotentials were applied to describe the electronic ex- 
change and correlation energy. The reaction path was described by the minimum energy path 
(MEP) between two local minimum states (Fig. Q]), and the MEPs for the desorption and diffusion 
of a hydrogen monomer on graphene were calculated using the climbing-image nudged elastic 
band method. 52 For calculating the vibrational frequencies of the periodic supercells, only the 
gamma point at the Brillouin zone center is selected. The charge-doped systems were compen- 
sated with the same numbers of opposite background charges. The optimized C-C bond length of 
graphene is 1.426 A; thus, a charge density of 1.0 x 10 13 cm" 2 corresponds to a charge number 
of 0.1314 e added into the graphene supercell. More details about the computational method and 
tests can be found in Ref. [281. 



III. RESULTS AND DISCUSSION 



Schematic drawings of the MEPs for the desorption and diffusion processes of a hydrogen 
monomer on graphene are shown in Fig. [TJ where the structures of the initial (reactant), SP (tran- 
sition) and final (product) states are also shown. The initial state is the chemisorption state, and 
the final state here is set as the physisorption state. The potential barrier is the energy difference 
between the SP state and the initial state. The variation of the adsorption energy of the hydrogen 
monomer in the initial state (E 1 ,), diffusion SP state (Ejjf) and desorption SP state (E d a e d f), and 
of the potential barriers for the desorption (AV des ) and diffusion (AVf^) of the hydrogen monomer 
with respect to the doping-charge density (<x) are shown in Fig. [21 

E l ads increases with the number of doping electrons (negative charge) or doping holes (positive 
charge). However, the curve of E 1 . is asymmetric about cr = 0.0, because the rate of increase 
is greater under hole doping than under electron doping. E d 'Jf increases monotonically with in- 
creasing cr from negative to positive, and the rate of increase becomes larger than that of E 1 , 
under hole doping, but smaller under electron doping. The curve of E d ^' s is somewhat flat com- 
pared with those of E 1 . and E'ff. The potential barrier of diffusion AV d '^ (= E ads - E dl Jf) 
decreases monotonically with increasing cr from negative to positive, while the potential barrier of 
desorption AV des (= E ads - E d d s f) increases with the number of doping electrons or holes. These 
phenomena in energy are related to the electronic structure of graphene and to the electron transfer 
and coulomb interaction between the hydrogen monomer and graphene, which will be explained 
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in detail below. 

The band structure of neutral graphene is shown in Fig. [3ta), where the anti-bonding n* band 
is above the Fermi level and the bonding n band is below the Fermi level. The anti-bonding and 
bonding characteristics of the n* and n bands are shown by their orbital shapes in real space in Fig. 
|3£b) and (c), respectively. Upon electron doping, the Fermi level goes up, leading to the occupation 
of the anti-bonding n* band, and the strength of the total IT bond is weakened. Upon hole doping, 
the Fermi level goes down, reducing the electron-occupation number in the bonding n band, and 
thus the strength of the II bond is also weakened. The E n obtained at different <x by Equ. @]are 
shown in Fig. Oa). E n decreases with increasing the number of doping electrons or holes, inverse 
to the variation of E 1 , . This is consistent with the weakening of the IT bond predicted from the 
electronic structure of graphene. The curve of En is asymmetric because the rate of decrease is 
greater under hole doping than under electron doping, like the variation of E 1 ^ . From Equ. |2]and 
|3j the weakening of the II bond, namely the decrease of E n , results in the increase of E ads and 
E* ds , both of which exhibit a similar trend (Fig. Eta)). However, the general lower value of E ads 
than that of E* ds indicates the effect of the charge doping on Ec-h, which is related to the electron 
transfer and coulomb interaction between the hydrogen monomer and the doped graphene (shown 
below). 

The calculated Lowdin electronic populations, also called the net atomic population by 
Mulliken, 53 of the hydrogen monomer (Nf) in the initial, diffusion SP and desorption SP states 
are shown in Fig. HI Nj? decreases monotonically with increasing cr from negative to positive. 
The Nj/s for the initial state are smaller than 1 .0 due to the electron transfer from the hydrogen 
monomer to graphene, mainly through the orbital hybridization between them. This orbital hy- 
bridization in the initial state will be shown later to be unchanged by the charge doping. However, 
the electron transfer is still influenced by charge doping due to the chemical-potential difference 
between the hybridized hydrogen monomer and graphene. The mono tonic decrease of Nf for the 
initial state indicates that the electron doping tends to block but the hole doping to promote the 
electron transfer from the adsorbed hydrogen monomer to graphene. This blockage of electron 
transfer under electron doping tends to reduce the affinity of graphene to bond with the hydrogen 
monomer, so that E C _ H decreases with increasing doping electrons, making E{. < E*. under 
electron doping (by Equ. © and ©). On the other hand, hole doping tends to promote electron 
transfer so as to increase E C _ H , making E 1 , a little larger than E*. at cr < 2.5 x 10 13 cm' 2 . How- 
ever, at cr > 2.5 x 10 13 cm" 2 , the further promoted electron transfer increases the coulomb repulsion 
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between the positively charged hydrogen monomer and the hole-doped graphene, which weakens 
the strength of the C-H bond and makes E 1 , < E*, again. The weakening of the strength of the 
C-H bond by charge doping can also be qualitatively observed from the variation of its bond length 
(cI C -h) under charge doping in Fig. [5J The variation of the protrusion height of the H-bonded C 
atom (Z c .) with a is also shown in Fig. [5J which will be used later. It should be noted that the 
increase of E 1 , by charge doping is mainly determined by the redution of the breaking energy 
of a total II bond (En), as described in detail in the previous paragraph, and the electron transfer 
and coulomb repulsion play minor roles. However, from the electronic density of states (DOS) 
of the neutral desorption and diffusion SP states as described in Ref. 28, the hydrogen monomer 
hybridizes little with graphene. This is the same for these doped cases, but only with Fermi levels 
tuned by charge doping (not shown). Thus, E de d s f and E d a 'Jf should be not significantly related 
to En- The Aff s for the desorption SP state are larger than 1.0 under electron doping and smaller 
than 1.0 under hole doping. The doping charges distribute across the system according to the re- 
quirement of electrostatic equilibrium, which explains the flatness of the E de d s f (a) curve in Fig. 
2fa). The N"s for the diffusion SP state are significantly smaller than 1.0. From the DOS in Ref. 
28L the hydrogen monomer will dope graphene with some number of itinerant electrons in this 



state. The monotonic decrease of (promotion of the electron transfer) with increasing <x is due 
to the chemical-potential difference between the hydrogen monomer and doped graphene, which 
is the same as the electronic population of the lithium adatom on charged carbon nanotubes. 54 
This promotion of electron transfer results in the increase of the affinity of graphene to bond with 
the hydrogen monomer in the diffusion SP state and results in the increase of E d 'Jf . As a result, 
the increase of AV des (= E ads - E de d s f) under charge (electron or hole) doping and the monotonic 
decrease of AV dlf (= E' ads - E d 'Jf) with increasing cr are conjunctly due to the effects of charge 
doping on the strength of the total II bond of graphene, the electron transfer and to the coulomb 
interaction between the hydrogen monomer and graphene. 

The calculated vibrational zero-point energy corrections (AF vib (0)) for the desorption and dif- 
fusion of H and D monomers under various <xs are shown in Fig. |6] The isotope effect of AF v $(0) 
is obvious; it decreases with increasing monomer mass. This will result in the reversed isotope 
effect in the activation energy by Equ. |6] AF vih (0) and its isotopic difference for desorption are 
generally larger than those for diffusion. These vibrational properties are related to the spectra of 
the localized vibrational modes of H and D monomers, which has been discussed in detail in Ref. 
28l The AF vjb (0)s for the diffusion of H and D monomers decrease with increasing cr, the same 
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as AVp'f. This is because the decrease of the potential barrier AVf^ will result in the decrease of 
the diffusion-MEP potential curvature around the initial state, which then results in the lowering 
of the frequencies of the effective vibrations for the diffusion of H and D monomers whose zero- 
point energies equal their A_F V! j(0)s. Except for the kinks at <x of -5.0 x 10 13 cm" 2 , the AF V! -&(0)s 
for the desorption of H and D monomers present the same variation as AV d p es with respect to cr, 
which is that the values increase with the number of doping electrons or holes. This is related to 
the increase of the potential barrier of the desorption-MEP and the increase of the frequencies of 
effective vibrations for the desorption of H and D monomers. The kinks in the AF v #,(0) curves 
for the desorption of H and D monomers at cr of -5.0 x 10 13 cm -2 are related to the structure of 
the desorption SP state. The protrusion height of the H-bonded C* (Z c .) in the chemisorption and 
SP states can be used to determine the structural information. It can be seen in Fig. [51 there is 
an obvious kink in the Z c .(cr) curve for the desorption SP state at cr = -5.0 x 10" 13 cm 2 , while 
the Z c .(cr) curves for the initial and diffusion SP states are much smoother. This kink makes the 
C* closer to the hydrogen monomer in the desorption SP state at cr = -5.0 x 10" 13 cm" 2 than 
at other charge densities. Then, the interaction between hydrogen monomer and graphene in the 
desorption SP state at cr = -5.0 x 10" 13 cm" 2 will be larger than at other charge densities, which 
also can be reflected in the spectra of the localized vibrational modes of the hydrogen monomer 
(shown below). 

The localized vibrational modes have large displacements (e(o>,), i is the index of the vibrational 
mode) of the hydrogen monomer in their eigenvectors of the vibrational dynamic matrix. The 
spectra of |e(a) ( )l 2 for the initial and desorption SP states at crs of -5.0 x 10 13 and -7.5 x 10 13 
cm -2 are shown in Fig. U\ The stretching (S) modes become imaginary in the desorption SP states. 
Compared with the bending (B) modes in the initial states, the bending modes in the desorption 
SP states are much less lowered at cr of -5.0 x 10 13 cm" 2 than at cr of -7.5 x 10 13 cm" 2 . The 
phonon spectra at other crs (e.g. see the Fig. 3 in Ref. 
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for the neutral case) are close to those 



at cr of -7.5 x 10 13 cm -2 and shift smoothly with cr, which is also reflected by the AF v a(0) curves 
(Fig. [6]). As described in the previous paragraph, the interaction between the hydrogen monomer 
and graphene in the desorption SP state at cr of -5.0 x 10 13 cm" 2 is larger than those at other crs, 
and thus (the absolute values of) the vibrational frequencies of the localized vibrational modes of 
the hydrogen monomer are higher at this cr than those at other crs. By Equ. [6l the result is that the 
AF V ib(0)s for the desorption SP state are lower at cr of -5.0 x 10 13 cm" 2 than at other crs, and that 
the E ac s for desorption at crs of -5.0 and -7.5 x 10 13 cm" 2 are very close to each other. By Equ. 
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|7]and from the detailed analysis in Ref. |28i the stiffening of the localized vibrational modes in the 
desorption SP state at cr of -5.0 x 10 13 cm" 2 will also result in the decrease of the corresponding 
v* m s compared with those at other <xs. 

The calculated jump frequencies (v) for the desorption and diffusion of H and D monomers on 
doped graphene are shown in Fig. [SJ The value of v decreases with increasing monomer mass. 
This isotope effect is due to the isotope effect in AF vi i,(0) (Fig. [6]) by Equ. [5]and[6l which has 



been discussed in detail in Ref. 
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In Fig. [8£a) and (b), the desorption v generally decreases with 
increasing the number of doping electrons (negative cr) or holes (positive cr), with the exception 
that the v for the desorption of H (D) monomer at cr of -5.0 x 10 13 cm -2 is less than that at cr of 
-7.5 x 10 13 cm" 2 . By Equ. [51 v is exponentially determined by E ac {cr) (= AV d p es + AF vjb (0)), where 
the linear dependence on v* qm usually plays a minor role. Except at cr of -5.0 x 10 13 cm" 2 , the 
E ac for the desorption of H (D) monomer increases with the number of doping electrons or holes; 
thus, the desorption v generally decreases with increasing the number of doping electrons or holes. 
However, the lowering of the AF v ^(0)s at cr of -5.0 x 10 13 cm" 2 makes the E ac s very close to those 
at cr of -7.5 x 10 13 cm -2 , and the relative magnitudes of the vs at these two charge densities are 
determined by the values of their v* m s. The v* qm at cr of -5.0 x 10 13 cm" 2 is smaller than those 
at other crs due to the stiffening of the localized modes (not shown), as described in the previous 
paragraph. In Fig. Hfc), the diffusion v increases monotonically with cr, which is due to the 
monotonic decrease of E ac for diffusion with increasing cr. Thus, it can be concluded that any kind 
of charge (electron or hole) doping will make the bonding between the hydrogen monomer and 
graphene kinetically more stable, and electron doping will prevent but the hole doping will trigger 
the diffusion of the hydrogen monomer on graphene. Thus, when increasing cr from negative to 
positive, there should be a crossover between the priorities of desorption and diffusion of H (D) 
monomer on graphene, like the crossover between AV p es and AVf^ in Fig. [2] 

In the annealing process, the system is kept at an annealing temperature (To) for a fixed time 
interval (f ). In the simulation here of the kinetic properties of H and D monomers in the anneal- 
ing process, to was set to be 300s. For electron-doped graphene layers, only the desorption of the 
hydrogen monomer was considered in the simulation, because the diffusion is stopped by electron 
doping. However, for hole-doped graphene layers, both the desorption and the diffusion of the 
hydrogen monomer were considered. The properties of desorption and diffusion can be character- 
ized by the relative residual monomer number (or the desorption probability) (n(to) /n(0)) and the 
diffusion radius (r d ,f(t )), respectively. The calculated variations of the n(t )/n(0)s and r di f(t )s of 
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H and D monomers at various crs with respect to Tq are shown in Fig. |9j At n(to)/n(0) = 0.5, 
the corresponding Tq was defined to be the desorption temperature. The desorption temperature 
of H (D) on neutral graphene is 282 (301) K. Under electron doping, the desorption temperatures 
of H (D) monomer are 295 (314), 333 (349) and 330 (348) K at crs of -2.5, -5.0 and -7.5 x 10 13 
cm -2 , respectively. Under hole doping, the desorption temperatures are 313 (331), 359 (376) and 
407 (424) K at crs of 2.5, 5.0 and 7.5 x 10 13 cm" 2 , respectively. The desorption temperatures for 
D monomer are about 18 K higher than those for H monomer. The r di f of H (D) monomer on 
neutral graphene at the desorption temperature is 0.6 (1.0) A, which is smaller than the C-C bond 
length of 1 .426 A and indicates the immobility of the hydrogen monomer in diffusion. Under hole 
doping, the r di fS of H (D) monomer at the desorption temperatures are 9.0 (12.5) nm, 0.65 (0.71) 
/im and 10.7 (9.8) /im at crs of 2.5, 5.0 and 7.5 x 10 13 cm" 2 , respectively. The isotope effect on the 
desorption and diffusion of the hydrogen monomer can be concluded to be that the lighter hydro- 
gen monomer is desorbed and diffuses more easily. The variations of the desorption temperatures 
and r di fS with respect to cr are the same as those of the vs for desorption and diffusion, respectively. 
If diffusion with radius above 0.1 //m is defined to be macroscopic diffusion, the macroscopic dif- 
fusion of H (D) monomer at temperatures below the desorption temperature can be achieved when 
cr reaches 5.0 x 10 13 cm -2 . 

In the constant-rate heating process, the system is heated at a constant rate (T = at). The 
heating rate (a) was taken to be 1 .0 K/s in the simulation here, which was a commonly used value 
in experiments. The calculated variations of the n(t )/n(0)s and r d! /(? )s at various crs are shown in 
Fig. [T0l The desorption temperature of H (D) on neutral graphene is 313 (333) K. Under electron 
doping, the desorption temperatures of H (D) monomer are 326 (346), 367 (384) and 363 (382) 
K at crs of -2.5, -5.0 and -7.5 x 10 13 cm" 2 , respectively. Under hole doping, the desorption 
temperatures are 345 (363), 394 (412) and 445 (461) K at crs of 2.5, 5.0 and 7.5 x 10 13 cm" 2 , 
respectively. The desorption temperatures in the constant-rate heating process are about 33 K 
higher than those in the annealing process. The desorption temperatures for D monomer are also 
about 18 K higher than those for H monomer. The r di f of H (D) monomer on neutral graphene 
at the desorption temperature is 0.8 (1.3) A, which is also smaller than the C-C bond length and 
indicates the immobility of the hydrogen monomer in diffusion. Under hole doping, the r di fS of H 
(D) monomer at the desorption temperatures are 7.5 (9.6) nm, 0.41 (0.41) /urn and 5.9 (5.6) fim at 
crs of 2.5, 5.0 and 7.5 x 10 13 cm" 2 , respectively. The isotope effect on the desorption and diffusion 
of the hydrogen monomer on graphene and the variations of the desorption temperatures and the 
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rjffS with respect to cr are the same as those in the annealing process. The macroscopic diffusion 
of the hydrogen monomer is also achieved in this constant-rate heating process when cr reaches 
5.0 x 10 13 cm 2 . 

Although the diffusion radius can be increased by hole doping, the hydrogen monomer on 
heavily hole-doped graphene layers may not be easily observed in experiments, because highly 
diffusive hydrogen monomers tend to meet and form hydrogen dimers and clusters, which are 
much more stable than isolated monomers, 3 or quickly diffuse to the edge of the finite-sized sam- 
ple. Thus, to observe the diffusion of a hydrogen monomer experimentally, a medium hole density 
was suggested. The hydrogen monomer on suspended graphene sheet has been detected to be 
stable by transmission electron microscopy (TEM) at room temperature. 21 Although this stability 
may be due to the low actual sample temperature caused by the cold trap used in the experiment, 
the effect of the incident electrons in TEM should also be considered. The incident electrons will 
inevitably excite the electrons in the bonding n band up to the anti-bonding n* band in graphene. 
This excitation reduces the electronic occupation number of the n band and increases that of the 
n* band together, which weakens the total n bond in graphene and enhances the stability of the 
hydrogen monomer on graphene, according to the analysis above. 

The adsorption of hydrogen monomer onto neutral graphene can break an aromatic n bond, 
which results in two dangling C(p z ) orbitals. Then, one C(p z ) orbital bonds with the H(ls) orbital, 
and the C(p z ) orbital left unsaturated forms an occupied spin-polarized quasilocal state (spin-up) 
around the hydrogen monomer . 23128,55 In the spectrum of the electronic DOS(Fig. [TTj) . the occupied 
quasilocal state presents as a narrow spin-up peak within the gap between the conduction and 
valence bands, and there is another narrow peak above the Fermi energy (E F ) that corresponds 
to the unoccupied spin-down quasilocal state. The energy difference between these two peaks is 
due to the exchange interaction of electrons, and this exchange splitting in energy was defined as 
A s . The electronic occupation number of these two quasilocal states can be changed by charge 
doping. The calculated DOS of the chemisorption state at various crs are shown in Fig. \TT\ The 
contribution of the H(ls) orbital to the total DOS (not shown) does not significantly vary with 
charge doping, and is nearly the same as that in the neutral case in Ref. 2M, which indicates 
that charge doping does not significantly influence the orbital hybridization between the hydrogen 
monomer and graphene. This can also be validated from the invariance of the gap of 1.29 eV 
between the conduction and valence bands with respect to cr. It can be seen that A s decreases with 
increasing \cr\. The electronic occupation numbers of the lower spin-up and the higher spin-down 
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quasilocal states were defined as n + and n~, respectively. It can be seen that under electron doping 
(cr < 0.0), the spin-up and spin-down quasilocal states are fully and partially occupied (n + = 1.0, 
rT < 1.0), respectively, while under hole doping (cr > 0.0), they are partially and not occupied 
(n + < 1.0, n~ = 0.0), respectively. The variations of A v and An (= n + - n~) with respect to cr are 
shown in Fig. [T3 The curve of A s (cr) is strictly linear and symmetric around cr = 0.0. A simple 
exchange model can be used to express A. s as 

A, = e - e + = -J; ff nS S " - {-J + eff )n + S + S + (15) 

where ± represent spin-up and spin-down quasilocal states, respectively; e ± are the energies of the 
quasilocal states; Jf f f are the exchange constants; 5* = ±j. In the exchange model, electrons 
with opposite spins do not interact with each other. The linearity and symmetry of the curve of 
A s (cr) indicate that the exchange constant does not vary with the spin orientation or cr, which can 
be expressed as -/j}y(cr) = Jeff( cr ) - Jeff - Then, Equ. [[5] can be rewritten as 

A * = \ J eff^ n ( 16 ) 

The magnetic moment of the system equals An fi B - The curve A i (cr) can be well fitted by Equ. [16] 
with J e ff taken to be 1.23 eV, as shown in Fig. [121 

IV. CONCLUSIONS 

The thermodynamic and kinetic properties of H and D monomers on doped graphene layers 
were studied using a composite method consisting of density functional theory, density functional 
perturbation theory and harmonic transition state theory. The calculated results were analyzed 
with reference to the electronic structures and phonon spectra. Electron doping has been found 
to heighten the diffusion potential barrier, while hole doping lowers it. However, both kinds of 
dopings heighten the desorption potential barrier. These phenomena in energies have been found 
to be conjunctly due to the effects of charge doping on the strength of the IT bond (defined to 
include the n and n* bonds) in graphene and to the electron transfer and the coulomb interaction 
between the hydrogen monomer and graphene. It has been found that hole doping is necessary 
for the observation of the diffusion of H and D monomers on graphene. The kinetic properties 
of H and D monomers on doped graphene layers during both the annealing process (to = 300 s) 
and the constant-rate heating process (a =1.0 K/s) were simulated. H monomer is more mobile 
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than D monomer in the kinetic simulations. Generally, both electron doping and hole doping can 
increase the desorption temperatures of hydrogen monomers. However, the diffusion of hydrogen 
monomers is prevented by the electron doping and triggered by the hole doping, and the diffusion 
radius increases with <x under hole doping. It has been found that the macroscopic diffusion of 
hydrogen monomers can be achieved at temperatures below the desorption temperature when cr 
reaches 5.0 x 10 13 cm' 2 in graphene. The effect of charge doping on the magnetic properties of the 
hydrogenated graphene were also studied. The exchange splitting of the spin-up and spin-down 
quasilocal states and the magnetic moment decrease linearly with the number of doped electrons 
(holes). The variation of the exchange splitting with respect to cr has been explained by a simple 
exchange model, where the exchange constant has been found not to vary with the spin orientation 
or the doping charge density. 
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FIG. 1: (Color online) The MEPs for the (a) desorption and (b) diffusion of the hydrogen monomer on 
graphene. The initial, saddle-point (SP) and final states in the paths are labeled with stars, and their struc- 
tures are shown alongside the graph. The yellow spheres are carbon atoms and the smaller blue spheres are 
hydrogen atoms. The final state in the desorption MEP is set to be the physisorption state. 




FIG. 2: (Color online) (a) The variations of E 1 ^ , E d ^f , E d 'Jf , En and E* ds with respect to cr. (b) The 
variations of AV des and Avf^ with respect to cr. 
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FIG. 3: (Color online) (a) The band structure of neutral graphene. The orbital shapes in real space of (b) 
the bonding n band and (c) the anti-bonding n* band at T, M and K points. In (a), the bonding n band and 
the anti-bonding n* band at T, M and K points are each labeled with a red circle. 




0.6 I 1 1 1 1 1 1 

-7.5 -5.0 -2.5 0.0 2.5 5.0 7.5 

a(10 13 cm" 2 ) 



FIG. 4: (Color online) The variations of N^s in the initial, desorption SP and diffusion SP states with respect 
to cr. 
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FIG. 5: (Color online) The variations of dc-H hi the initial state and Zc*s in the initial, desorption SP and 
diffusion SP states with respect to cr. 
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FIG. 7: (Color online) The calculated spectra of |e(w;)| 2 (i is the index of the vibrational mode) of H (D) 
monomer in (a) the initial state and (b) the desorption SP state at cr of -5.0x 10 13 cm -2 , and in (c) the initial 
state and (d) the desorption SP state at cr of -7.5 X 10 13 cm -2 . 
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FIG. 8: (Color online) The variations of the jump frequencies for the (a, b) desorption and (c) diffusion of 
H and D monomers on graphene at various crs with respect to the inverse of temperature. 
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FIG. 9: (Color online) The variations of the n(t Q )/n(0)s and r^,/(7())s for H and D monomers at various crs 
with respect to the annealing temperature (To) in the annealing process (to = 300s). The r c nf(to)s at negative 
crs are not shown, because the diffusion of H (D) monomer is stopped by electron doping. 
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FIG. 10: (Color online) The variations of the n(t)/n(0)s and rdif(t)s for H and D monomers at various crs 
with respect to temperature in the constant-rate heating process (a = 1.0 K/s). The r^if(t)% at negative crs 
are not shown, because the diffusion of H (D) monomer is stopped by electron doping. 
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FIG. 1 1 : (Color online) The electronic DOS for the chemisorption states at various o~s. The gray vertical 
dotted line at the Fermi level guides the eyes. The definition of A s is shown in the DOS spectrum at cr = 0.0. 
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FIG. 12: (Color online) The variation of An, A v and jJ e ffAn with respect to cr. J e ff is fitted to be 1.23 eV. 
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